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We have studied the competition between helix formation and aggregation for a simple polymer 
model. We present simulation results for a system of two such polymers, examining the potential 
of mean force, the balance between inter and intramolecular interactions, and the promotion or 
disruption of secondary structure brought on by the proximity of the two molecules. In particular, we 
demonstrate that proximity between two such molecules can stabilize secondary structure. However, 
for this model, observed secondary structure is not stable enough to prevent collapse of the system 
. into an unstructured globule. 
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Aggregation of partially folded proteins is a problem in many bio-related fields. In bioprocessing operations, 
mechanical and/or environmental stresses lead to protein unfolding which may result in irreversible aggregation and 
loss of protein activity?! Another problem in bioprocessing is the recovery of partially-folded recombinant proteins, 
produced in densely packed "inclusion bodies". Recovery of these proteins requires dissolving the bodies in denaturant 
and then slowly diluting the solution, minimizing unwanted aggregation such that the protein folds to the native 
state. Finally, protein aggregation has been implicated in the origin of amyloid fibrils, which have been linked to 
many different diseases? 2 ^" 4 - Understanding the process of aggregation, and determining the conditions under which 
it will take place, is therefore a matter of some importance. 

The competition between protein folding and protein aggregation is believed to be determined by the balance of 
inter and intramolecular forces. Protein folds are stabilized by strong intramolecular interactions between hydrophobic 
and hydrogen bonding groups in the protein interior. Aggregation of folded proteins is of necessity weak, mediated 
by interactions between heavily solvated charged and polar surface groups. However, at high temperatures, or upon 
' addition of denaturant, the fold becomes less stable. This structural disruption causes hydrogen bonding and hy- 
. drophobic groups to become exposed, providing "sticky" sites for aggregation. The resulting competition between 
\Q) ' structure formation and aggregation may hold the key to understanding the process of aggregation. 

The details of this process are difficult to experimentally resolve, since partially folded intermediates have very 
low solubilities. As such, computational approaches are extremely useful. Detailed, all-atom models have been used 
\ to study the process of aggregation.^! 7 -^ However, such models are computationally intensive, and details of the 
force field used can cloud underlying behavior. Simple "reduced models" that can produce protein-like secondary 
| structur o 10 ' 11 ' 12 ' 13 ! 14 ! 15 ! 16 ! 17 ' 18 are therefore useful as a parallel approach to the problem. 

In a previous paper^ we proposed an isotropic polymer model which can form helical secondary structure. This 
model consists of a chain of overlapping spherical beads, and is related to the well-studied tangent-sphere square- well 
polymer chain model. Low temperature transitions to "crystalline" ordered states have been observed for tangent 
sphere polymer o 20 i 21 i 22 Further, interactions between pairs of tangent-sphere polymers are well studied? 2 ^ 

In this paper, we extend the study of the overlapping bead model by examining the behavior of isolated pairs of 
such polymers. In particular, we examine how the formation of helical structure affects the interactions between the 
two molecules and influences the structure of the resulting complex. The remainder of this paper is organized as 
follows. In section [ill we describe the simulation methodology and the particulars of the model. In section HlH we 
present the results from the simulations, and in section HVt we present our conclusions. 
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c5 ' II. SIMULATION DETAILS 



The polymer model we consider consists of a linear chain of N impenetrable spheres of diameter a. The distance 
between bonded spheres is restricted to vary between 0.9Z and l.ll, where I is the nominal bond length. Directly 
bonded spheres may overlap (i.e., a > I), although non-bonded spheres may not. In addition to excluded volume 
interactions, the spheres also interact via a square-well potential. The interaction potential u(r) between the spheres 
is given by: 

{oo for r < a 
— e for a < r < Act (1) 
for Act < r 
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Table I: Temperatures used in simulations. 
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where r is the distance between the centers of the spheres, e is the strength of the square-well interaction, and A 
characterizes the width of the square-well interaction. In this work, A = 1.5. Since this potential is pairwise, the 
interaction energy of a system of polymers can be decomposed into an intramolecular energy Ei n t ra (interactions 
between monomers on the same molecule) and an intermolecular energy Ei nter (interactions between monomers on 
different molecules). At ajl — 1.0, the model reduces to the traditional tangent sphere polymer model; interactions 
between pairs of such polymers has been studied extensively elsewhere i 23 i 24 i 25 i 26 ' 27 For larger values of a/l, we have 
observed helix formation™, despite the isotropic, achiral nature of the model. While the precise mechanism of helix 
formation for this model remains unclear, it has been suggested that helices are in fact the optimal packing structure 
for any string-like object^ in the same fashion that optimal packing of spheres yields hexagonal or cubic close-packed 
periodic lattices. 

Monte Carlo simulations were performed for a pair of N — 20 (20mer) square- well molecules with a/l = 1.3, 1.6, and 
1.9. Pivot, crank-shaft, and sphere displacement moves were used to sample the individual polymer conformations. 
Overall translation and rotation moves of the polymers were performed to sample the relative positions and orientations 
of the polymers. Umbrella sampling of the polymer separation was performed with a potential of the form 

U(r) = a{r - bf (2) 

where r is the center of mass separation between the polymers, and the parameter a = 40fcsT was used to sample 
intervals in polymer center of mass separation. The parameter b was varied from a value of to 25?. The values of b 
were initially distributed evenly in units of I, and extra values were added wherever insufficient overlap in r histograms 
was observed. For each umbrella potential, parallel tempering in temperature was used to enhance equilibration of 
the low temperature systems. The temperatures used in the simulations are given in Table U . For each umbrella 
potential, the multiple histogram method was used to combine simulation data at different temperatures. 

Each production simulation run consisted of a total of 10 6 parallel tempering (box swap) moves. Between each of 
these attempted moves, 80 pivot, crank-shaft, and sphere displacement moves were attempted, as well as 80 molecular 
translation and rotation moves. Samples of the system configuration were taken after every 10 3 attempted parallel 
tempering moves. For each set of conditions, an equilibration run of 10 5 attempted parallel tempering moves was 
performed, followed by 10 production runs. The reported results are averages over the production runs, and the 
uncertainties were estimated by the jackknife method. 

We are particularly interested in the influence of intermolecular interactions on the helical structure which has been 
observed for this model in isolation^ As such, structures are characterized through the torsion r. For a continuous 
space curve R(s), the torsion is defined as: 



(R'(s)xR"( s ))-R"'(s) 

t(s) = 5 (3) 

|R'(s) x R"(s)| 2 

where primes denote derivatives with respect to s. For the polymer chain, the required derivatives are calculated 
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using a central difference approximation: 

R" = Ri+i — 2Ri + Rj_i 

R-i" = 2 (R*+2 - 2 (Ri+i - Rj_i) - R»-2) (4) 

where Ri is the vector position of monomer i. An instantaneous molecular conformation is characterized by the 
configurational average torsion f , given by: 

N-2 

f =Y J n/{N-A) (5) 

i=3 

where Tj is the torsion at monomer i, defined on i £ [3, N — 2], and the bar denotes an instantaneous configurational 
average. Chiral (torsional) symmetry breaking is used to detect helix structure formation; helix formation is associated 
with a probability distribution p (f ) which is bimodal, symmetric, and with p (f = 0) w 0. Combined with observation 
of helical snapshot configurations, this provides strong evidence for helicity. Structures are also characterized by the 
absolute of the mode of p (f), denoted f max . 

It is also instructive to consider statistical correlations between properties of the molecules. The correlation between 
observables A and B is given by: 



m (A m {{AB) {A) {B))2 (*\ 

cor (A, B) = —. r— r- (6) 

(y?) ^f) (V) - <£> 2 ) 

where angle brackets denote ensemble averaging. A correlation of zero indicates that the two variables are com- 
pletely uncorrelated, while a correlation of one indicates perfect linear correlation. We measure correlations between 

the intramolecular energies of the individual molecules cor (^E^ tra , E^ tra ^j , the torsions of the individual molecules 

cor(fi,T2), and the orientations of the end-to-end bond vectors, (|? e (i) ' ^e(2)|)- The last quantity should take a value 
of one if the end-to-end vectors of the molecules are aligned, zero if the end-to-end vectors are perpendicular, and a 
half if the orientations are uncorrelated. The absolute value of the dot product is taken because the molecules have no 
directionality; antiparallel orientations are equivalent to parallel orientations. For unstructured globules, we expect 
this quantity to be one half (since the end-to-end vector should be randomly oriented), whereas for helical structures, 
the end-to-end vectors are a reasonable approximation for the helical axes, and the dot product provides a measure 
of the orientational correlation of the helices. 

The pair potential of mean force (PMF) ip(r) is estimated from the observed histograms h(r) of separation between 
molecules. For a given simulation, this is given by: 

ip{r) - -k B Tln(h(r)/r 2 ) - U(r) +C (7) 

where fee is the Boltzmann constant and C is an unknown energy offset, arising from the incomplete normalization 
of probability histograms from simulation. The U(r) term removes the bias introduced by the umbrella potentials. 
Since solvent is implicit, i^(r) — at separations for which the polymers can have no direct interaction (i.e.. r/l > 
(N — 1) + \a/l). For simulations sampling these large values of r, it is possible to fix C by this criterion. Results 
for neighboring umbrellas are then "patched" together by finding values for C which minimize the total statistical 
uncertainty in the combined results. 

The last quantity of interest is the second virial coefficient B 2 , which may be calculated from the PMF via: 



B 2 = 2tt J (1 - cxp (-/3ip(r))) r 2 dr (8) 

Positive values for B 2 indicate an overall repulsive force acting between the molecules; negative values indicate an 
overall attractive force. 
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Figure 1: Single molecule phase diagram, 
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Figure 2: Contour plots of the estimated probability distributions p (r; r) for a/l = 1.6 at (a) ksT/e — 0.255, (b) fesT/e — 0.435, 
(c) k B T/e = 0.651, and (d) k B T/e = 1.164. 



III. RESULTS AND DISCUSSION 



In previous work^ we have studied this polymer model for isolated 20mers, down to a temperature of fcsT/e = 0.2. 
It was found that, at sufficiently low temperatures, these polymers may form helical structures, depending upon the 
value for a/l. Further, two separate helical phases were observed over the range 1.475 < a/l < 1.66. For reference, 
the phase diagram of the single polymer system is reproduced in Fig. [TJ In this work, we study the behavior of an 
isolated pair of polymer molecules with a/l = 1.3, 1.6 and 1.9. 



A. a/l = 1.6 

We begin with a description of the results for the a/l — 1.6 system, which best demonstrates behaviors which 
also occur in the a/l = 1.3 and a/l = 1.9 systems. Results are presented across a range of separations at four 
representative temperatures; ksT/e = 0.255, for which isolated polymers are in the helix 2 phase, ksT/e = 0.435, for 
which isolated polymers are in the helix 1 phase, ksT/e = 0.651, just above the helix 1-globule transition temperature, 
and ksT/e = 1.164, at which the polymer is well into the globule phase. 

In order to illustrate the structural behavior of the interacting molecules, contour plots of p (f; r) against r are 
displayed in Fig. [2 The thermodynamic behavior is shown in Fig. [3l The potential of mean force is given in Fig. [3fa) , 
the mean intramolecular interaction energy (Ei ntra } is given in Fig. [3](b) , and the variation of the intermolecular 
interaction energy (E inter ) is given in Fig. ^c). The correlations between the intramolecular energies cor (E\, E 2 ), 
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Figure 3: Variation of thermodynamic properties with molecular separation for cr/l = 1.6: (a) potential of mean force, (b) 
intramolecular interaction energy, and (c) intermolecular interaction energy. Solid lines denote ksT/e = 0.255, dotted lines 
denote ksT/e = 0.435, dashed lines denote ksT/e = 0.651, and dot-dashed lines denote ksT/e = 1.164. Inset: Root mean 
square radius of gyration r a /l for fcsT/e = 0.255. Errors are less than symbol size. Equivalent data at higher temperatures are 
indistinguishable at this scale. 
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Figure 4: Intermolecular correlations for a/l = 1.6 between: (a) intramolecular energies, (b) conformational torsions, and 
(c) end-to-end bond vectors. Solid lines denote fcsT/e = 0.255, dotted lines denote fcsT/e = 0.435, dashed lines denote 
ksT/e — 0.651, and dot-dashed lines denote fcsT/e = 1.164. 



the mean torsions cor (n, f 2 ), and the orientation of the end-to-end distance vectors (r e m • r e(2)), for the molecules 
are shown in Fig. [H 

The lowest representative temperature (fc^T/e = 0.255) is less than the helix 1-helix 2 transition temperature, and 
the results are typical of this region. In Fig.[2ja), it can be seen that at separations r/l > 7.5, p(f;r) is independent of 
r, bimodal and approximately zero at f = 0. In Fig. [3] (solid line), it can be seen that, at these large separations, ip(r) 
and {Ei n ter) are approximately zero, and (Ei ntr a) appears r-independent. Finally, as shown in Fig. [4] (solid line), all 
three measured correlators are approximately zero for large separations. This behavior indicates that, for r/l > 7.5, 
the two molecules do not significantly interact with each other, existing independently in the helix 2 state. There is 
a weak interaction, since observed (Ei nter ) is non-zero, with magnitude greater than observed error, for r/l < 8.5; 
however, these interactions do not appear strong enough to significantly perturb the isolated equilibrium state. 

The interactions between the molecules begin to perturb their individual structures at r/l rs 7.5. In the range 
4 % r/l < 7.5, p(f;r) remains bimodal, however, the positions of the peaks, ±? max depend upon separation. As 
shown in Fig. [31 over the same range, ip(r) and (E int er) are monotonically increasing functions of r, while (E intr a) 
remains approximately constant. Further, the torsions and end-to-end vectors of the molecules become strongly 
correlated in this region, while intramolecular energies remain uncorrelated (see Fig. [4j) . 

The bimodality of the torsion distributions, as well as the approximately constant value for the intramolecular 
energies, indicates that the molecules maintain helical configurations within this range. A snapshot configuration of 
two helical molecules taken at this temperature is presented in Fig. UK a) . The observed correlations show that the 
helices have a strong preference to take on the same handedness, and lie end-to-end along the same axis. 
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(a) (b) 

Figure 5: Snapshot configuration for (a) a/l = 1.6, ksT/e = 0.255, umbrella potential centered at b/l = 5, and (b) a/l = 1.6, 
fcsT/e = 0.651, umbrella potential centered at b/l — 5.6. 

The average end-to-end distance for the unperturbed polymer at this temperature 19 is r e /l ~ 5; two completely 
rigid "average" helices would only begin to interact at r/l = 5.5. For two helices to interact at r/l = 7.5, they must 
have r e /l > 7, which corresponds to an extension of at least 27%. The interactions between the molecules appear to 
stretch the molecules, relative to their isolated equilibria. 

The stretching effect weakens upon closer approach between the molecules; at separations 4 < r/l < 5, f max 
approximately returns to its value at large separation. Further, (|r e (i) •r e ( 2 )|) decreases to approximately 0.5 at 
separations in this range, indicating that the two molecules are less likely to align. For an isolated molecule at this 
temperature, r e /l w 5; for two aligned "average" helices, there would be a longitudinal compression at separations 
closer than this. As such, the stretching force is eliminated, and in order to remain helical, the helical axes must 
lose some orientational correlation upon further approach. Note that the torsion correlator remains high across this 
region. In the range 4 < r/l < 7.5, the intramolecular energies of the two molecules are uncorrelated. 

Once the two molecules approach closer than r/l « 4, p (f ; r) becomes unimodal (see Fig. [2]); that is, the structures 
become achiral, and helicity is lost. This behavior is marked by a small discontinuity in total energy (data not shown) 
at r/l — 4.15 and a sudden loss of torsional correlations between the molecules. For 2.5 < r/l < 4, i/j(r) and {Ei nte r) 
increase monotonically with r, (Ei ntra ) remains approximately constant, and no significant correlation exists between 
the intramolecular energies. The intermolecular interactions in this range are strong enough to fully disrupt the helical 
structures, causing both molecules to adopt individual unstructured globule configurations. We suggest that the lack 
of intramolecular energy correlations and the approximately constant intramolecular energy indicate that, at these 
separations, the complex consists of two globules with only limited interpenetration. The steady decrease in (Ei nte r) 
with decreasing r suggests that the two globules deform upon approach to give a larger contacting interface. 

At r/l — 2.5, however, there is a large discontinuity in {E int er) and in {E int ra)- Further, the slope of (E int ra) versus 
separation changes at this separation; for r/l > 2.5, the average intramolecular energy is approximately constant 
(with respect to r), however, for r/l < 2.5, it is a decreasing function of r. For r < 2.5, there is a negative correlation 
between the intramolecular energies. This shift has no obvious signature in p(f;r). The inset to Fig. [3] shows that, 
below r/l k 2.5, the radius of gyration of the system becomes approximately constant (r g /l m 2.1); further approach 
does not lead to contraction of the complex. The total energy of the system (i.e., E — Ei ntra + Ei nter ) also becomes 
approximately constant with respect to r at these short separations (data not shown). This suggests increasing 
interpenetration between the two molecules, with intramolecular contacts being replaced on a one-to-one basis by 
intermolecular contacts. Since ip(r) is also approximately constant for r/l < 2.5, this implies that the entropy of the 
system also remains approximately constant at these separations. Strong anticorrelation of intramolecular energies, 
indicated by large negative values for corr {E\, E2), shows that any "swelling" (loss of intramolecular contacts) of one 
molecule must be matched by an equivalent "collapse" (gain in intramolecular contacts) of the other. These data are 
consistent with two interpenetrated molecules in a collapsed, unstructured globule state. 

It is unclear why, for separations 2.5 < r/l < 4, the two molecules act like immiscible droplets, whereas for r/l < 2.5, 
the molecules can interpenetrate; nor is it clear why the transition between the two is so sharp. 

Qualitatively similar behavior is seen when ksT/e — 0.435, which is in (and typical of) the helix 1 phase for an 
isolated polymer. Once again, at large separations, molecules adopt two independent, uncorrelated helical structures. 
Upon closer approach, interactions once again cause correlation between the helical axis and torsions and cause 
stretching forces which lead to extended helical configurations. Closer still, correlations and helical structure are 
suddenly lost, and the system appears to form two contacting but distinct unstructured globules. Finally, interactions 
cause the two molecules to interpenetrate, forming a single globule, which again is marked by a broad, flat minimum 
in ip(r). 

Increasing the temperature to ksT/e = 0.651, there is unexpected behavior. This temperature is above that 
of the helix 1-globule transition, and, consequently, at large separations the molecules adopt unstructured globule 
configurations, with no chiral symmetry breaking. However, at r/l — 6.8, there is a small discontinuity in the total 
energy (not shown), and p(f;r) changes from a broad, unimodal distribution to a bimodal distribution. That is, 
upon approach, interactions between the two molecules cause chiral symmetry breaking and helix formation, at a 
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Figure 6: Contour plot of the absolute mode Tmax of p (r; r, T) for a/l = 1.6. The positions of the globule - helix 1 and helix 1 
- helix 2 transitions are marked on the right hand side. Note that nOIl _ Z6rO Tmax is not sufficient for chiral symmetry breaking; 
see discussion in Section |n] and Ref. 19. 

temperature above that where helical configurations would normally be stable. A snapshot configuration is shown in 
Fig. [U[b) . The character of this pair of helices is very similar to the helices seen for the same separations at lower 
temperatures (see Fig. [U[a)), with similar torsion values and correlation in helical axes and torsions. We suggest 
that attractive intermolecular interactions between the globules produce a stretching force, which has previously been 
shown to stabilize helix formation.-^ As the two molecules continue to approach each other, there is once again a loss 
of helical structure and intermolecular correlations to two unmixed globules, followed by interpenetration and mixing. 

Finally, consider the temperature ksT /e = 1.164, which is above that of the helix-globule transition for isolated 
polymers. The torsion distribution p (f; r) exhibits a broad, unimodal structure for all separations; this system shows 
no chiral symmetry breaking, and no helix formation. There are no observed correlations between the torsions or the 
end-to-end vectors of the molecules. There is an apparent transition between separate and interpenetrated globules, 
marked by a slight discontinuity in total energy, changes in the slope of ip(r) and the inter- and intramolecular energies, 
and sudden onset of strong anticorrelations between intra-molecular energies. It is also interesting to note that the 
potential of mean force at this temperature is qualitatively similar to those for the lower temperatures; it appears 
that the presence or absence of helical structure has little effect on the character of the PMF. 

A contour plot of f max as a function of /csT/e and r is shown in Fig. [6l Note that regions of non-zero f max do 
not necessarily indicate helical configurations; as pointed out in our previous work-i& for single molecules, there exists 
a region just above the helix transition where p(f) is double peaked but has a non-zero value around f = 0, and as 
such is not truly symmetry breaking. However, this contour plot does schematically show the increase in stability of 
helical structures at intermediate separations, as well as the loss of helical structure at close separation. 



Potential, torsion and correlation data for the a/l = 1.9 system are shown in Figs. El [8] and [9l A plot of f max 
against separation and temperature is shown in Fig. [TUl Once again, we consider representative temperatures: this 
time, ksT/e = 0.651, typical of the region just above the helix 2 - globule and ksT/e = 1.455, typical of higher 
temperatures. We were not able to fully equilibrate the two molecules at temperatures below the helix 2 - globule 
transition. 

The behavior at the lower temperature feT/e = 0.651 is qualitatively similar to the situation for the a/l = 1.6 
system at /csT/e = 0.651. The torsion histograms in Fig. [UJa) reveal a range of separations over which helical 
structures are observed at temperatures higher than found for isolated molecules. This range is associated with a 
correlation between the torsions of the two molecules (Fig. [9](b) , solid line) , although this correlation is not as strong 
as that observed for the a/l — 1.6 system. Interestingly, there is no significant correlation in the end-to-end bond 
vectors over this range (Fig. EDJc), solid line). At close separations, the intramolecular energies of the two molecules 
become strongly anticorrelated (Fig. [9](a) , solid line) . This correlation appears gradually on approach from the loss 
of helical structure (at r/l « 4) to a discontinuity in intramolecular and intermolecular energies (at r/l = 2.4). Once 
again, this separation coincides with the point at which the radius of gyration (not shown) becomes approximately 
constant, with r g /l w 2.4, and below this separation, total energy and ip(r) (Fig. [71(a), solid line) are also independent 
of r. Aside from the lack of orientational correlation, the only other qualitative difference is a repulsive "hump" in 
ip(r) at long ranges for the a/l = 1.9 system. The maximum of this barrier (r/l — 8.25) is near to the separation 
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Figure 7: Variation of thermodynamic properties with molecular separation for a/l = 1.9: (a) potential of mean force, (b) 
intramolecular interaction energy, and (c) intermolecular interaction energy. Solid lines denote ksT/e = 0.651 and dotted lines 
denote k B T/e = 1.455. 
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Figure 8: Contour plots of the estimated probability distributions p(r;r) for a/l = 1.9 at (a) ksT/e = 0.651, and (b) 
fcsT/e = 1.455. 



at which (Ei nte r) becomes non-zero — this barrier is, therefore, entropic in nature, corresponding to loss of entropy 
of extended configurations. The long range of this hump suggests that such extended configurations are relatively 
frequent occurrences — this is to be expected, considering the stiffness of the a/l — 1.9 polymer. Overall, then, the 
results are consistent with the interpretations offered for the a/l — 1.6 system at the same temperature. 

At ksT/e = 1.455, there is new behavior. While the torsional behavior (Fig. [8Kb)), potential of mean force and 
energies (Fig. [TJ dotted line), and the torsional and end-to-end vector correlators (Fig. [9](b) (c), dotted line) appear 
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Figure 9: Intermolecular correlations for a/l — 1.9 between: (a) intramolecular energies, (b) conformational torsions, and (c) 
end-to-end bond vectors. Solid lines denote fcsT/e = 0.651, dotted lines denote ksT/e — 1.455. 
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Figure 10: Contour plot of the absolute mode T ma x of p(r;r) for a/l = 1.9. All temperatures shown are above the helix 2 - 
globule transition temperature. 

similar to those seen at high temperatures for a/l = 1.6, there are strong internal energy correlations across the range 
2 r S 8 (Fig. Efa), dotted line). These correlations are observed at temperatures up to ksT/e = 2.273 - above 
this, the system returns to the qualitative behavior seen for a/l = 1.6, fcsT/e = 1.164. Snapshot configurations taken 
across this range of temperature and separation do not show any remarkable common features; the system appears 
to be exploring a large ensemble of unstructured configurations. 

The positive intramolecular energy correlation indicates that swelling of one polymer is associated with swelling 
of the other. We relate this to the stretching force that we have inferred to operate between helices over a similar 
range of separations. If one of the polymers is in a particularly open (high intramolecular energy) configuration, 
then it will have a larger radius. This will increase the likelihood of interactions with the other polymer; if these 
interactions are attractive, then the other polymer can only expand (at fixed intermolecular separation). The variance 
in intramolecular energies for the a/l = 1.9 system in the globule phase is approximately double that observed for the 
a/l = 1.6 globule phase (data not shown). This indicates that the a/l = 1.9 globule phase has much lower internal 
cohesion than the globules at smaller overlaps, which we suggest is due the inherent stiffness of the polymer, and 
which allows cross interactions to produce these cross-correlations which are not observed for the a/l = 1.6 system. 

In summary, the behavior of the a/l = 1.9 system is largely similar to that of the a/l = 1.6 system, with destruction 
of helical structure at short range, stabilization of helical structure at intermediate range, and an apparent transition 
between separate and interpenetrated globules at very close range. A summary of the torsional behavior of the system 
is shown in Fig.[10l Two significant differences in behavior occur. First, strong correlations in intramolecular energies 
between the two molecules are observed in the unstructured globule phase. We suggest that this does not occur in the 
a/l = 1.6 system because the individual globules are more internally cohesive, as indicated by lower intramolecular 
energy fluctuations. Further, no correlations between the end-to-end bond vectors of the helices are observed; we will 
return to this point in the conclusion. 



Finally, consider square-well polymers with a/l = 1.3. Here, we present two representative temperatures, fcsT/e = 
0.255 and ksT/e = 0.651. Single 20mers with a/l = 1.3 do not display*! 9 any helical structure down to a temperature 
of ksT/e = 0.2. Potential of mean force and energies are shown in Fig. [TTJ contour plots of p (f ; r) in Fig. HH and 
correlations in Fig. [T3l 

At the lower temperature ksT/e = 0.255, the behavior is similar to that seen for the a/l — 1.6 and a/l = 1.9 systems 
at temperatures just above the helix-globule transition. At long range, there is no torsional symmetry breaking. This 
continues on approach down to r/l = 5.65, where there is a sudden discontinuity in intramolecular and intermolecular 
energies. At this separation, p(f; r) abruptly becomes bimodal, and the torsions and end-to-end vectors become 
strongly correlated. This behavior indicates that the helix phase is stabilized by interactions between the polymers. 
This stability continues across the range 3.65 < r/l < 5.65. At r/l — 3.6, there is a small discontinuity in the total 
energy, the torsion histogram is no longer bimodal, and the torsions and orientations of the molecules are no longer 
correlated. At closer separations, there is a gradual increase in anticorrelation between the intramolecular energies 
up to a discontinuity in intermolecular energy at r/l = 1.8. Upon closer approach, the intramolecular energies are 
strongly anticorrelated, and the PMF and total system r g remain approximately constant. This behavior is once again 
interpreted as a transition from two immiscible polymer globules at intermediate separation to two interpenetrating 
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Figure 11: Variation of thermodynamic properties with intermolecular separation for a/l = 1.3: (a) potential of mean force, 
(b) intramolecular interaction energy, and (c) intermolecular interaction energy. Solid lines denote ksT/e = 0.255, and dotted 
lines denote fcgT/e = 0.651. 
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Figure 12: Contour plots of the estimated probability distributions p(r;r) for a/l = 1.3 at (a) ksT/e — 0.255, and (b) 
k B T/e = 0.651. 



globules at close separations. 

At the higher temperature ksT/t = 0.651, the behavior is similar to that seen for the a/l = 1.6 and a/l = 
1.9 systems in the unstructured globule regime. The torsions and end-to-end vectors remain uncorrelated at all 
separations, and the torsion histograms remain unimodal. The mixing transition between the two globules appears 
at close separation, marked by the onset of intramolecular energy correlations. No strong intramolecular energy 
correlations are observed. The variance in the intramolecular energies is of the order of magnitude of that seen for 
the a/l = 1.6 system, which suggests that the globules have enough internal cohesion to prevent such correlation. 

The torsional data are summarized in Fig.[l4]and are similar to the data presented in Figs. [6] and [10] for temperatures 
above the helix-globule transition. 



D. Summary 



In summary, we have found similar behaviors for pairs of model polymers for a/l = 1.3, 1.6 and 1.9. In each 
case, intermediate separations between the molecules enhance helix stability, showing paired helical configurations 
at temperatures above the single molecule helix-globule transition temperature. For these separations, the paired 
helices have strongly correlated chiralities, and for a/l = 1.3 and 1.6, are aligned end-to-end. At separations closer 
than the typical helix length, however, the molecules lose helicity and adopt disordered configurations. At separations 
immediately below the typical helix length, the disordered globular configurations of the polymers are uncorrelated. As 
the separation becomes approximately equal to the radius of gyration of the complex, however, energy anticorrelations 
begin to develop between the two molecules; below this separation, the radius of gyration of the complex remains 
approximately constant. The potential of mean force is qualitatively insensitive to a/l, consisting of a monotonically 
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Figure 13: Intermolecular correlations for a/l = 1.3 between: (a) intramolecular energies, (b) conformational torsions, and (c) 
end-to-end bond vectors. Solid lines denote fcsT/e = 0.255; dotted lines denote fcsT/e = 0.651. 
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Figure 14: Contour plot of the absolute mode T m ax of p(r;r) for a/l = 1.6. 



decreasing "ramp" from the initial onset of interactions at large separations, down to the separation at which energy 
correlations begin to occur between globules. Below this separation, the PMF remains approximately constant in a 
broad, flat minimum to zero separation. For the temperatures at which helical configurations are observed, the PMF's 
are all strongly attractive, as evidenced by very large negative values for the second virial coefficient £?2, shown in 
Fig.Ql 

We suggest that the stabilization of helical structure at moderate separation occurs because attractive interactions 
between the two molecules exert a stretching force. Stretching forces have been shown to stabilize helix formation 
in other simple helix- forming models^ Correlation in chirality and helical axis can be seen as a consequence of 
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energy minimization. By sampling face-to-face configurations, contacts between helices are possible at intermediate 
separation; by having the same chiralities, the match-up between the faces, and hence the number of contacts, is 
maximized. For a /I = 1.9, there is no significant axial alignment, but there is a strong correlation in chirality between 
helices. This suggests that face-to-face interactions are strong enough to promote identical chirality between molecules 
even when such interactions are not strong enough to hold the two polymers in axial alignment. 

IV. CONCLUSIONS 

The observed behaviors are similar to those seen in simulations of helical homopolypeptides^ up to the separation 
at which helical structure is lost. Upon close approach, simulated homopolypeptides are observed to switch to parallel 
paired helical configurations, whereas for the simple model presented here, the system prefers to collapse into an 
unstructured globule. The helical structures which the system exhibits are simply not stable enough to overcome the 
free energy to be gained by formation of the collapsed structure. 

One may therefore ask how the helix can be stabilized against the collapsed, disordered globule. We consider that 
the 0-point of a polymer (the temperature at which B2 = 0) is related to the gas-liquid transition in a bulk system. 
The 9-point is also the temperature below which a single polymer will begin to collapse from an extended, disordered 
coil state (analogous to a gas) to a collapsed, disordered globule state (analogous to a liquid). If we hypothesize that 
the helical configuration is analogous to a bulk crystalline phase for the polymer model (since it is an ordered state 
which must be reached through a first-order transition) , then we can consider methods by which the crystalline phase 
may be stabilized against the liquid phase. It has been show n 30 ' 31 that, for bulk systems of particles interacting via 
a simple, isotropic potential, decreasing the range of the attractive part of the potential can decrease the stability of 
the liquid phase to the point of leaving the liquid-gas critical point metastable with respect to the crystal phase. In 
analogy, we suggest that by decreasing A in Eq. ([I]) , we may depress the 0-point of the system without depressing the 
helix transition temperature, so that the stability of the dense globule may decrease to the point that the observed 
collapse may not occur, giving results closer to those from simulations of homopolypeptides. 

This suggestion may have a basis in the chemistry of real homopolypeptides. The strongest interaction present for 
these molecules is likely to be hydrogen bonding between segments of the peptide backbone. This hydrogen bonding 
is very strongly anisotropic, occurring for only a small range of orientational alignments between peptide segments. It 
has been shown that for bulk systems interacting via simple anisotropic potentials, decreasing the range of orientations 
over which interactions can occur depresses the temperature of the liquid- gas critical point.— 

We also observe that moderate proximity between pairs of square well homopolymers can promote the stability of 
helical structures. We have suggested that this may be linked to the stabilization of helical structure upon stretching 
of a single molecule, observed for other simple helix-forming models.— For pairs of molecules, this stabilization is not 
enough to prevent collapse into an unstructured globule. Experimentally, helix forming proteins aggregate to form 
ordered structures with molecules in aligned, extended configurations. Since such configurations are not seen for this 
model nor for simulations using pairs of more detailed model peptides,- it seems that they must be stabilized by three 
or more body interactions between molecules. This suggests that bulk simulations using the model presented here 
may be useful. 

In conclusion, we have shown that interactions between pairs of square well homopolymers can promote helix 
stability, while at close proximity, interactions cause collapse from paired helices to a single disordered globule. We 
have further suggested a modification to the model (a reduction in the range of the attractive part of the potential) 
which may stabilize the helix phase against such collapse, in order to provide better qualitative agreement with 
simulations of pairs of helix forming homopolypeptides. Finally, we have suggested that the reason for the formation 
of amyloid aggregate structures may lie in many-body interactions between molecules. These results underline the 
need for a better understanding of the underlying nature of the helix-disordered globule transition, in order to gain 
an understanding of how interactions between molecules can increase or decrease the relative stabilities of the various 
structures. 
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